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We consider two classes of helical magnets. The first one has magnetic ordering close to antifer- 
romagnet and the second one has magnetic ordering close to ferromagnet. The first case is relevant 
to cuprate superconductors and the second case is realized in FeSrOs and FeCaOa. We derive the 
effective field theories for these cases and calculate corresponding excitation spectra. We demon- 
strate that the "hourglass" spin-wave dispersion observed experimentally in cuprates is a fingerprint 
of the "antiferromagnetic spin spiral state" . We also show that quantum fluctuations are important 
for the "ferromagnetic spin spiral" , they influence qualitative features of the spin- wave dispersion. 
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I. INTRODUCTION 

Helical spin systems have being studied for long time. 
Physics of these systems is fascinating and sometimes 
is very much different from that of collinear spin sys- 
tems. The noncollincar spin systems manifest quite un- 
usual temperature-driven phase transitions and critical 
properties,^"'* for a review see Ref. 5. Excitation spectra 
(magnons) in noncoUinear systems are also very different 
from that in collinear magnets, see, e.g., Ref. 6. Generi- 
cally there arc two ways to generate a noncollincar spin 
structure, spin frustration and Dzyaloshinsky-Moriya 
spin orbit interaction. In the present paper we consider 
only the spin-frustration mechanism. It is worth to note 
that the Dzyaloshinsky-Moriya mechanism leads to in- 
teresting quantum field theories,'''® but this is outside of 
the scope of the present work. 

The usual field theory approach to a spin helix induced 
by spin frustration is based on the following representa- 
tion of the spin S, see Ref. 5 

S{t, r) = a{t, r) cos(Q • r) + b{t, r) sin(Q • r) . (1) 

Here Q is the wave vector of the spin spiral and a{t,v), 
b{t, r) are the dynamic fields. The nonlinear Ginzburg- 
Landau action describing temperature-induced phase 
transitions is written in terms of the fields a and b, see 
review.^ Typically a frustration implies that Q ^ 27r/a, 
where a is the lattice spacing. Typical momenta, that 
appear in fluctuations of the a and b fields, arc much 
smaller than Q, and hence Q is not influenced by the 
fluctuations. 

Sometimes, due to accidental fine tuning of parame- 
ters, the wave vector of the spin spiral is small, Q <^ 
27r/a, i.e. the system is close to a collinear ferromag- 
net. This happens in FeSrOa and FeCaOa compounds.^ 
In this case typical momenta of fluctuations are compa- 
rable with Q, and hence the standard approach^ based 
on the Ginzburg-Landau action for a and b is not valid. 
Another example is the case of Q close to the vector 



of antiferromagnetic ordering. On the two-dimensional 
(2D) square lattice this means Q (7r,7r) -I- Q. This 
is the "antiferromagnetic spin spiral" realised in cuprate 
superconductors. ^°'*^ In this case Q is small because of 
the small doping x, and hence, once more, the approach 
based on Eq.(l) is not applicable. 

There is one more important difference of the field the- 
ories considered in the present paper from the case (1). 
This difference concerns a direction of Q. In the cases 
we consider, the vector Q, in spite of being small, is al- 
ways correlated with lattice directions. For example for a 
square lattice it can be directed only along a diagonal or 
along a side of the square. This implies an essential cor- 
relation between the short range dynamics. Ax ^ a, and 
the long range dynamics. Ax ^ Hence, this also 

implies that 0(3) and S0(3) nonlinear a-models cannot 
be applied for these situations. 

In the present paper we consider effective quantum 
field theories for helical spin magnets close to antiferro- 
magnets and to ferromagnets. The structure of the paper 
is the following. Section II reminds the reader the logic 
of derivations of the effective action for collinear antifer- 
romagnets and ferromagnets. In Section HI we consider 
the "antiferromagnetic spin spiral" and demonstrate that 
the "hourglass" spin- wave dispersion observed in neutron 
scattering from cuprates is a fingerprint of such state. In 
Section IV we consider the "ferromagnetic spin spiral" 
and show that quantum fluctuations influence qualitative 
features of the spin-wave dispersion. Section V presents 
our conclusions. 



II. COLLINEAR ANTIFERROMAGNETS AND 
FERROMAGNETS 



To derive the effective action for a collinear isotropic 
antiferromagnet (AF) one can use the following hcuris- 
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tic arguments.^ (i) The staggered magnetization is de- 
scribed by a unit vector n, fp = 1. (ii) The elastic energy 
is proportional to (Vn)^. (iii) The Larmor theorem: If 
there is a state of the system without magnetic field, then 
the same state in a uniform magnetic field B must pre- 
cess with frequency lo — Bp' (iv) The static energy can 
depend on the magnetic field only quadratically because 
the system has no net magnetization. In other words, n 
is not a vector, it is an axis. The only Lagrangian that 
satisfies the above requirements is the Lagrangian of the 
nonlinear cr-model 



C 



AF 



{Vnf 



(2) 



n — n X i? 

2 V ^ V 2 

where the coefficients are the perpendicular magnetic sus- 
ceptibility x± E^nd the spin stiffness ps ■ 

To derive the effective action for a coUincar ferromag- 
net one can use similar arguments, (i) The magnetization 
is described by a unit vector n, = 1. (ii) The elastic 
energy is proportional to (Vn)^. (iii) The Larmor theo- 
rem: If there is a state of the system without magnetic 
field, then the same state in a uniform magnetic field B 
must precess with frequency lo = B. (iv) The static en- 
ergy depends on the magnetic field linearly because n is 
true magnetization. The only Lagrangian that satisfies 
the above requirements reads 



^ = 77^(777 . 



SiB-0 



(3) 

Here S is spin of the ferromagnet and 77 is a CP^ two- 
component spinor. The field ^ is a bosonic field, we just 
use the CP^ representation for the field. Note that the 
difference between antiferro and ferro is only in the point 
(iv): quadratic or linear dependence of energy on an ex- 
ternal uniform magnetic field. 

It is worth reminding that Lagrangians and fields in 
a quantum filed theory always depend on normalization 
point. Throughout the paper we use normalization at 
zero momentum (infinite distance). This implies that in 
the groimd state the normalization is |(77)| = |(^)| = 1. 
This also implies that spin stiffnesses and the magnetic 
susceptibility in (2) and (3) are renormalizcd parameters 
that include all quantum fiuctuations. These are not bare 
parameters. 

In subsequent sections we derive and analyse effective 
actions that describe spin spiral states close to antifer- 
romagnetic and ferromagnetic cases. This means that 
the wave vector of the corresponding spiral is small com- 
pared to the inverse lattice spacing. Specifically we have 
in mind cuprates for the AF case and FeCaOa for the 
ferro case. 



^ A coUinear state implies that the minimum dimensionality of the 

system is 2D + time. 
^ Hereafter we set g/ig = 1, g is the g-factor, and fiB is Bohr 

magneton. 



III. SPIN SPIRAL CLOSE TO 
ANTIFERROMAGNET. BOSONIC FIELD 
THEORY INSPIRED BY CUPRATES 

A. Effective action 

Here we consider the 2D case (2D -I- time) , the 3D gen- 
eralization is straightforward. The effective action for the 
spin spiral state in the generalized t — J model has been 
derived in Ref. 11. The action is written in terms of 
the bosonic field ft, describing staggered spins of the par- 
ent Mott insulator, and the fermionic field ipa describing 
mobile holons. The spinor ipa describes pseudospin of 
the holon, the pseudospin originates from two AF sub- 
lattices. The index a = a,b (fiavor) describes a nodal 
pocket where the holon is located. The action reads 



C = 



(4) 



/3 



ipiD^ipa + V2g{ipiaiPa) ■ [n x (e„ • V)r 



The first line in the Lagrangian corresponds to the non- 
linear cr-model (2) describing spins. The second line rep- 
resents time derivatives of fermionic fields. The first term 
in the third line is the kinetic energy of fermions, /? is the 
inverse mass, the second term in the third line was de- 
rived long time ago by Shraiman and Siggia.^° The last 
line represents the interaction of the pseudospin a with 
an external uniform magnetic field. -'^^ The unit vectors 
ea = (I/V2, 1/V2) andefa = (1/^2, -l/\/2) are directed 
along the corresponding diagonal of the square lattice. 
The original t-J model is gauge invariant. Therefore the 
effective action (4) is also gauge invariant and this is re- 
flected in covariant derivatives 



Dt = dt + -a ■ [ft X n] 
D = V + - [nx Vn] 



(5) 



acting on the fermionic field. 

The Lagrangian (4) describes many physical effects: 
superconductivity, spin spiral, particle-hole decay, even 
phase separation for some values of parameters. In the 
present work we aim only at the spin spiral. Therefore we 
will truncate fermionic degrees of freedom in (4) throw- 
ing away all physical effects except of the spin spiral. We 
retain only spin and pseudospin degrees of freedom re- 
sponsible for the spin spiral. The way of truncation is 
pretty straightforward, we replace the fermionic field ipa 
by a bosonic field describing pseudospin in the CP^ 
representation 



(6) 
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Hence, the truncated Lagrangian reads 
Caps - ^ (k - [n -k B]f - {V nf 



(7) 



where is the effective "spin" of the ^-field, the stiffness 
of the ^-field is p^. Tlie coupling constant between n 
and fields is g, it docs not coincide with g in (4), but 
they are related. Importantly, the truncated Lagrangian 
(7) still respects the gauge invariance of the original t-J 
model since it contains only covariant derivatives acting 
on the pseudospin, Dt is defined in Eq.(5), and the spacial 
derivative is defined as 



(8) 



A related important point is that (7) respects also the 
Larmor theorem. If there is a state without magnetic 
field that satisfies equations of motion, then in a magnetic 
field this state precesses with frequency uj = B. This is 
a sort of a global Ward identity. 

It is worth noting that statuses of the Lagrangians 
(4) and (7) are different. The Lagrangian (4) has been 
derived from the generalized t-J model in a controlled 
way^^. The small parameter that controls accuracy of 
the derivation is small doping. On the other hand we 
obtain the Lagrangian (7) from (4) by an uncontrolled 
truncation of fermionic degrees of freedom. As a result 
we get an effective Ginzburg-Landau action (7) that con- 
tains only spin n and pseudospin ^. In the present paper 
we analyze only the action (7) which, as we show below, 
describes experimental observations remarkably well. 

The conventional approach^ to the spin spiral state is 
based on Eq. (1), it contains two vector dynamic fields 
a(t, r) and h{t, r). Our approach also contains two vector 
dynamic fields n{t,r) and £,{t,r). The difference is that 
our effective action (7) is valid for any g ^ 1, it is not 
restricted by the region q <^ Q. We will show that in 
the limit q <^ Q our approach results in two indepen- 
dent Goldstone modes with linear dispersions. Hence, 
as one should expect, in this limit both approaches are 
equivalent. 

Hereafter we set the external magnetic field equal to 
zero, B ~ 0. Neglecting quantum fiuctuations one can 
easily find the ground state of the Lagrangian (7), 

n = (C,5,0) 

C ~ cos(Q • r) 

S = sin(Q-r) 

fa = 6 = (0,0,1) 

g ea + Bb 



Q 



V2 



-(1,0) 



(9) 



Here we choose x- and y-directions in the spin space to lie 
in the plane of the spin spiral. Both and £,b are directed 
along z in the spin space, and in this case the wave vector 
of the spiral Q is directed along the (1,0) direction of the 
square lattice. One can also choose = —^b = (0, 0, 1). 
In this case Q oc — eb is directed along the (0, 1) 
direction of the square lattice. The choice between these 
two possibilities is spontaneous. 



B. The in-plane spin- wave excitation 

Consider the spin-wave excitations above the ground 
state (9). There is the in-plane excitation (p = ip{t,r), 
where 

n = (cos(Q • r + (p), sin(Q ■ r + ip), 0) 

Ca- (0,0,1), (10) 

and the out-of-plane excitation n^, where 



n = (C, S,0)^/l-nl + (0,0, n^) 



(11) 



Note that we use indexes x, y, z to denote axes in the spin 
space and we use indexes 1,2 to denote directions in the 
coordinate space, the direction 1 is parallel to the spin 
spiral vector Q and the direction 2 is perpendicular to 
the spin spiral vector. 

We first consider the in-plane magnon (10) that docs 
not involve fiuctuations of the ^-ficld. Substituting ex- 
pression (10) in the Lagrangian (7) and keeping only 
quadratic terms in the small perturbation ip one finds 



X-L .2 



(12) 



Hence the Green's function of the in-plane magnon is the 
same as the Green's function of a magnon in the parent 
CT-model, 



oj^ — c^q^ 



iO 



where 




(13) 



(14) 



is the bare spin-wave speed. The in-planc magnons can 
be excited in inelastic neutron scattering. Neutrons in- 
teract only with the n-field that represents real spins, 
therefore the scattering amplitude is given by the Fourier 
transform of H from (10) and the scattering probability 
for unpolarized neutrons with the energy transfer w and 
the momentum transfer k is proportional to 



/i„(a;,k) 



1 

2^ 



Im [G^(w, k - Q) + G^(w, k + Q)](15) 
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1 Pi 



FIG. 1; The spin wave dispersion of the in-plane magnon as 
it is seen in neutron scattering 



Hereafter we introduce the dimensionless energy and the 
dimensionless momentum 



UJ 



P 



(16) 



The momentum is measured in units of the wave vector 
of the spin spiral and the energy is measured in units of 
the magnon energy in the parent cr-model with momen- 
tum of the magnon equal to the spin spiral wave vector. 
According to Eq.(15) the in-plane magnon is seen in neu- 
tron scattering as two cones originating from the points 
±Q. The corresponding plot for the direction along the 
spin spiral is shown in Fig. 1. 



C. The out-of-plane spin-wave excitation 

We have shown in the previous subsection that the in- 
plane magnon in the spin-spiral state remains trivial. It 
is not the case for the out-of-plane magnon. Substitution 
of (11) in Eq.(7) gives the following action expressed in 
terms of and ^^j^ 



(Vn.)^] (17) 

a ^ 

[(V^a. + Q^cy + SQn, + CVn^f 



Ei 

2 



1 



9 
V2 



We keep here only quadratic terms. Let us perform in 
(17) two subsequent gauge transformations. The first one 



Cax 



exp<i --(Q-r-)(T^ 



^^ay 



(18) 



The second transformation reads 



cxp 

Saa; 



Say ^Oiy 



(19) 



Geometrically the gauge transformations (18) and (19) 
describe the transition to the proper reference frame of 
the vibrating spin spiral. In terms of the new variables 
r/" and £" the Lagrangian (17) consists of three parts 



^AFS — -Cji + + Cint 

r X± ■ 2 Ps i-r-T \2 r^'2 2 



(20) 



Hereafter we skip the double prime in notations of ry" 
and The first part, £„, describes the free motion 
of the -field; the second part, describes the free 
motion of the ^-fields; and the third part, describes 
the interaction between and £. Equations of motion 
generated by the action (20) are 
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E 



V2 



(e„ • V)Cay+ pdQ-^Ko 



= -2 



Pi^iay 



PsQi' 



^ay 



-^{e^ ■ V)n, - p^Q ■ V)n, 



P^^Cax - '^"^ Co 



(21) 



Looking for solution in the form 



<,ax — J-*a^ 

t _ _ -iujt+iq-r 
<,ay — tv^ac 



(22) 



we find the relations between the coefficients A, Ba, Ca, 
and the frequencies of excitations. 



^Ip^ ^ip = liOi + O2) T \j\{0^ - O2Y + O3 

01 = p2 + 52 

02 = a'(l + foV)' 

2 2,2 

P =Pl+P2 ■ 



(23) 
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These equations are written in terms of dimensionless 
frequency and momentum (16). We remind that pi is 
parahel to Q and P2 is perpendicular to Q. In (23) we 
have introduced two dimensionless parameters 



QPs 
cS^ 




(24) 



Note that these parameters are unrelated to the a,b- 
flavors. 

Importantly, there are two branches of the dispersion 
that originate from two different degrees of freedom, spin 
ft and pseudospin ^. In order to reproduce the spin-wave 
excitation spectra of the t-J model at small doping we 
chose b << 1, a > b. The parameter b defined in (24) 
determines the top of the lower branch ili.p, and the pa- 
rameter a defined in (24) determines the bottom of the 
upper branch fl2.p- Plots of Qi^p and fl2.p, for directions 
of p parallel to the spin spiral and perpendicular to the 
spin spiral are presented in Figs. 2A,B for values of pa- 
rameters b = 0.25, a = 0.5. Similar plots for b ~ 0.25, 
a = 0.3 are presented in Figs. 3A,B. In Figs. 2 and 3 the 




FIG. 2: (Color online) A and B: the dispersion Qfc of the out- 
of- plane magnon in directions parallel to the spin spiral, (A); 
and perpendicular to the spin spiral, (B). Values of parame- 
ters of the effective field theory are a — 0.5, b — 0.25. The 
blue dotted lines show the magnon dispersion in the parent 
cr-model. Thickness of Qi and f22 lines is proportional to the 
corresponding quasiparticle residue Zi, Z2. The quasiparticle 
residues are plotted separately for the same directions in Figs. 
C and D. 

dispersion is shown by very thick lines. The thickness is 
not a decay width, the decay width in this approxima- 
tion (no alive fermions) is negligible. The thickness is 
proportional to the corresponding quasiparticle residue 
calculated below. 

To find absolute values of the coefficients A, Ba, Ca, 
in (22) we need to quantize (20). There are two magnon 
creation operators, cj^, and c\p, corresponding to two 
branches of the dispersion. The standard quantization 




FIG. 3: (Color online) A and B: the dispersion fife of the out- 
of-plane magnon in directions (A) parallel to the spin spiral, 
and (B) perpendicular to the spin spiral. Values of parameters 
of the effective field theory are a = 0.3, b — 0.25. The blue 
dotted lines show the magnon dispersion in the parent a- 
model. Thickness of Qi and Q2 lines is proportional to the 
corresponding quasiparticle residue Zi, Z2. The quasiparticle 
residues are plotted separately for the same directions in Figs. 
C and D. 



gives 

= ^ Ni^p + h.c. 
1=1,2 



(25) 



1 



i,p 



^ ani^p [V2{ea ■ p) + b'^pi] 
ic^x^- 2^ P5 Ni,p + h.c. 

. «^(l + &V)[V2(e„.p)+6V] 

^ay = I 2^ TT2 75 Ni^p + h.c. 

Here index ^ = 1,2 enumerates the dispersion brunches, 
the expression for O2 is given in (23), and Zi ^ is the 
quasiparticle residue 



z 7 (O1-Q2) 



.(26) 



The quasiparticle residues Zi p and Z2.p are plotted in 
Figs. 2 C,D and in Figs. 3 C,D. Interestingly, the quasi- 
particle residues are practically isotropic, the dependence 
on the direction of p is extremely weak. The magnon 
Green's function defined in the standard way reads 



G„Jw,9) = -i J {T {n,{r,t)n,{0,0)})e''^*~"'''dtd'^r 



(27) 



We remind that w, and q, p are related according to 
Eq.(16). The inelastic neutron scattering intensity for 
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unpolarizcd neutrons with the energy transfer to and the 
momentum transfer fc is proportional to 



/o„t(w,k) = -ilm G„^(w,k) 

TT 



(28) 



Figs. 2 and 3 show the "hourglass" dispersion 
that agrees with the typical dispersion observed in 



The physical meaning of relations (30) and (31) is very 
simple: the reduction of static response is transferred to 
the dynamic response. The expressions (30) and (31) 
must be renormalized by subtraction of the ultraviolet- 
divergent contribution that corresponds to the parent o- 
model. Since the Green's function of the in-plane exci- 
tation (13) is exactly the same as the magnon Green's 



cuprates. 13-15 The lower branch has strongly function in the parent cr-model, the fluctuation ((^2) 



anisotropic momentum dependence with zeros at the 
Goldstone points q — iQ, while the upper branch 5^2, q 
is practically isotropic in the momentum space quickly 
approaching the bare spin-wave dispersion at larger q. 
The intensity (quasiparticle residue) of the lower branch 
Z\^q is strongly peaked at q = 0, at the "neck" of the 
"hourglass" , and the intensity decays dramatically away 
from this point, see Figs. 2 C,D and Figs. 3 C,D. On 
the other hand the intensity of the upper branch Z^^q 
is zero exactly at q = and quickly approaches unity 
away from this point. The dispersion obtained in the 
present work demonstrates that the "hourglass" shape 
is a "fingerprint" of the spin spiral driven by the inter- 
action between spin and pseudospin. The two branches 
ill and ^1 reflect the two degrees of freedom: spin and 
pseudospin. 



D. Quantum fluctuations 

The spin spiral state that we consider can be char- 
acterized by the vector order parameters, that indicate 
ordering of spin n and pseudospin ^, and a pseudoscalar 
order parameter V that indicates chirality. The order 
parameters are 



(el) , (6) 



~V2 



V)n\) . 



(29) 



In semiclassical approximation values of the vector order 
parameters are given by Eqs. (9), the chiral pseudoscalar 
is 7-" = Q. Note that P changes sign under time reflection. 

Due to quantum fluctuations values of some order pa- 
rameters are reduced. The reduction of the spin field 
is 



|(n)| = (v/l-^2_„2)^l 



1, 2 

2\ z 



(30) 



Only the length is reduced, the spacial dependence is still 
given by Eq. (9), (n) = |(n)|(cos(Q • r),sin(Q • t),0). 
To calculate the reduction of the length we note that 
the expectation values (tp^) and (n^) can be expressed in 
terms of the Green's function. 



E 

q 



27ri 



Gm(a;,q) 



(31) 



renormalized to zero. The out-of plane fluctuation gives 



Zi.^ 



1 

2^ 



(2^)2 



(32) 



Here the subscript R stands for "renormalized" , X — > In 
. The dimensionless integral can be easily calculated, 
it depends only on the parameters h and a defined in 
(24). Values of the integral for a and h corresponding to 
Figs. 2,3 are 



J/{(a = 0.5,6 = 0.25) = 0.23 
Ji?(a = 0.3,6 = 0.25) = 0.15 



(33) 



The pseudospin is not a gauge invariant object, it 
cannot be directly measured and its expectation value de- 
pends on gauge. One can calculate the expectation value 
of defined in the original action (7) or one can calcu- 
late the expectation value of obtained after the gauge 
transformation (18), (19). Below we calculate ((^^.l)^)' 
but we skip the double prime in notations. Using Eqs. 
(25) we find 



J = 



E 

;=i,2 



Ps 



(34) 



62)2 



20. 



i.p 



At p — > oo the integrand in is approaching the con- 
stant value equal to a. This corresponds to the ultraviolet 
divergence of {^'^) in the parent ^-niodcl defined by Eq. 
(3). Therefore must be renormalized. — > JTr, in the 
following way 



Jr = I {(■••) 



4 



(35) 



where (...) stands for the integrand in (34). After the 
renormalization JTr is convergent. 

The chiral pseudoscalar order parameter V is defined in 
Eq. (29). In this definition the field ^q, is the pseudospin 
of the original action (7), and in this sense the chiral 
order parameter is gauge invariant. It is easy to rewrite 
the order parameter in terms of the gauge transformed 
^" and n^. 



V 



1 



1 



= l-7^{<)--{ii^)- 



1 



V2Q 



(36) 
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where we again omit the double prime in notations. The 
first two terms after the unity have been already cal- 
culated in Eqs.(32),(34), and the last term immediately 
follows ft-om (25), 



(37) 



E 

1 = 1.2 



Zi^p 2a2[l + 6V]r + <Pp 



20, 



The integrand in /C decays at large p as oc so /C is 
ultraviolet divergent. We will discuss this divergence in 
subsection E. 

Quantum fluctuations do not influence direction of the 
spin spiral wave vector Q, but they change the absolute 
value of the vector. To calculate the quantum fluctuation 
correction to Q we need to account for nonlinear terms in 
the Lagrangian (7). To do so it is convenient to generalize 
the notations (10) and (11) for the in-plane and out-of- 
plane excitations in the following way 



n 




h (0,0,1)5', 




~ i£,ax, Cain \J ^ ^ 






~ cos(Q ■ r + ip) 






= sin((5 • r + (ys) 




C-y 


= cos 7 




iS'-y 


= riz = sin 7 . 





(38) 

The gauge transformation (18), (19) is generalized as 
= cxp|-^(Q-r + (y5)cr2|exp|^7crj^|7?^ 

Substitution of this transformation in (7) gives the 
quadratic Lagrangian (20) as well as the following cubic 
term (again we omit the double prime in ^) 




(40) 



By calculating the tadpole with if as an external leg, see 
Fig. 4, we flnd 




a 

Value of the wave vector of the spin spiral follows from 
the condition that there are no linear in ViyS terms in the 
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FIG. 4: The tadpole diagram giving the quantum correction 
to the spin spiral wave vector. 



effective action. The corresponding term in (7) in the 
leading approximation is —ps{Q-V(p) + -^ J2ai^a-'^'f)- 
The condition of the disappearance of this term gives the 
value of Q presented in Eq. (9). The account of the 
tadpole (41) gives the quantum correction to the spiral 
vector, 

f = Q-^^)(-.!)4(fi.) 



a 

The last term here is very shuilar to (37), 

^(EWQ-v)..)^g/c' 



(42) 



(43) 



^ Zi,p 2a2[l + 6V][l + 6^]p? d^P 
^f^^ 2^;:; [ill^ - O2] (2^ 



It is also ultraviolet divergent. We discuss the divergence 
in the next subsection. 



E. Discussion and conclusion of the 
"antiferromagnetic spiral" section 

Both IC in Eq.(37) and /C' in Eq.(43) are uftravio- 
let divergent. There is no a formal way to renormalize 
out these ultraviolet divergences. The effective bosonic 
fleld theory deflned by the action (7) is unrenormalizable. 
Practically this means that quantum corrections depend 
on the ultraviolet cutoff and different physical quantities 
depend on this cutoff in different ways. So it is impos- 
sible to reduce renormalization of all physical quantities 
to the renormalization of a finite set of parameters. 

Thus, the simplest conclusion is that within the ef- 
fective low energy bosonic theory (7) one can calcu- 
late the spin-wave spectra. The results are presented 
in Figs. 1,2,3. On the other hand it is impossible to cal- 
culate quantum corrections. However, we can somewhat 
improve the situation if recall that the effective action (7) 
originates from the effective action (4) that is equivalent 
to the t-J model at low doping. Therefore, the ultra- 
violet cutoff for fluctuations of the n-field (32) is about 
the inverse lattice spacing, A„ ~ 1/a. This cutoff does 
not appear in our results explicitly as soon as we use the 
renormalized spin stiffness and spin-wave velocity of the 
Heisenberg model^^ ps = 0.18J, c — 1.65J (we set the 
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lattice spacing equal to unity, a = 1). The ultraviolet cut- 
off for fluctuations of the ^field (34), (37), (43) is of the 
order of the Fermi momentum of holes, ^ pp ^ ^/t^, 
where a; ^ 1 is the doping level. At q > pf descrip- 
tion of fermions in terms of a structureless bosonic field 
^ does not make sense. The wave vector of the spin spi- 
ral scales linearly with doping^-'^, Q (x x. Therefore, wc 
conclude from Eqs. (37) and (43) that the quantum fluc- 
tuation corrections to the spiral wave vector and to the 
pseudoscalar order parameter are 

^ oc 6^7^ , ^ oc +Vi . (44) 

The ^/x dependence originate from the integrals K, and 
/C' that are regularized by imposing the ultraviolet cut- 
off ~ ^/x. There are two important points about 
Eqs. (44). (i) The quantum correction to the spiral wave 
vector is proportional to b^. According to experimental 
spectra^^"^^ b < 0.25 and hence 6^ < 0.06. So, we pre- 
dict a very small quantum correction to the vector of the 
spin spiral, (ii) The sign of the SV is positive since /C 
in (37) is negative. Hence quantum fluctuations enhance 
the pseudoscalar order parameter. This trend is opposite 
to the trend in (n). According to Eq. (32) quantum fluc- 
tuations reduce value of |(ri)| and at a sufficiently large 
Q the static component of spin vanishes, |(ri)| = 0. Dis- 
appearance of static spin does not mean that the spin 
spiral disappears, the spin spiral just becomes dynamic. 
Hence, |(n)| = indicates the Quantum Critical Point 
(QCP) for transition from the static spin spiral to the 
dynamic one. The spin spiral itself is characterised by 
the pseudoscalar chiral order parameter V independently 
whether this is static or dynamic case. This order param- 
eter goes smoothly through the QCP. So, the spin spiral 
survives in the magnetically disordered phase. 

There arc two types of spin wave excitations, the in- 
plane magnons and the out-of-planc magnons. The spec- 
trum of the in-plane magnons is given by Eq.(15) and 
Fig. 1, and that of the out-of-plane magnons is given 
by Eq.(28) and Figs. 2,3. It follows from Eq.(15) that 
neutron scattering with excitation of in-plane magnons 
is necessarily accompanied by diffraction on the static 
spiral. Therefore, one should expect that the intensity 
of the inelastic neutron scattering with excitation of in- 
plane magnons is sharply decreasing at approaching the 
QCP. The intensity is zero in the magnetically disor- 
dered phase. On the other hand the intensity of neu- 
tron scattering with excitation of out-of-plane magnons 
(Figs. 2,3 with the hourglass dispersion) is not related 
to a diffraction and hence it is not very sensitive to 
the QCP. The two branches of the hourglass disper- 
sion originate from two different degrees of freedom, spin 
and pseudospin. The lower branch of the hourglass has 
strongly anisotropic momentum dependence with zeros 
at the Goldstonc points q — ±Q, while the upper branch 
is practically isotropic. This behavior agrees with the 
typical dispersion observed in cuprates.^'^^^^ 



IV. SPIN SPIRAL CLOSE TO FERROMAGNET. 
THE FIELD THEORY INSPIRED BY FeSrOa 

The present section is inspired by neutron scattering 
data^ from FcSrOa and FeCaOa that indicate the spin 
spiral close to a collinear ferromagnet. The wave length 
of the spin spiral in FeSrOs is about 10 lattice spacing. 
The compound is an insulator, therefore there is only a 
spin field in the game, there is no a pseudospin, and there 
is no other relevant dynamic variable. To understand 
the kinematic structure of the effective action we refer to 
the minimal Heisenberg model suggested in Ref.^^ In the 
present paper we consider generic scenarios, therefore we 
address for simplicity the 2D case. The 3D generalization 
is straightforward. 

A. Effective action 

The Hamiltonian of the minimal Heisenberg model is 
H — — Jl Si ■ Sj + J2 Si' ■ Sj' + J3 Si" ■ Sj" . 

ij i'j' 

(45) 

determined on the square lattice shown in Fig. 5. The 
Jl link is ferromagnetic and J2, J3 links are antiferro- 
magnctic. For our purpose it is sufficient to consider a 




FIG. 5: 2D Heisenberg Ji — J2 — Ji model on square lattice. 
There is spin S at each site. 

classical spin, 5 3> 1. Let us impose the spiral configu- 
ration 

5^= S'(cos(Q • r),sin(Q • r),0) . (46) 

The classical energy per unit cell corresponding to the 
Hamiltonian (45) reads 

E = ^jS*^ {- Ji[cos((52;a) + cos((5^a)] 

-t- 2 J2 cos((5:ca) cos(Qj,a) 

4- J3[cos(2Q^a) -t-cos(2Qj^a)]} , (47) 

where a is the lattice spacing. Notice that in the present 
section we denote by indexes x and y directions in the 
coordinate space. Assuming that Qa <^ 1 and expanding 
the energy in powers of Q we find 

E = -^Q-^{2J2+U3-Ji) (48) 
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(I6J3 + 2J2 - Ji)a^ 
24 



To have a small nonzero Q we need a small negative spin 
stiffness term. Thus, the effective energy per unit cell in 
terms of the ferromagnetic n-field, ~ 1, reads 



E{n) = -ips(Vn)2 



+ 6i(9,a,n)2 + -62[(aln)2 + (9,2,^)2] 



(49) 



The parameters are 



= ^2(2J2+4J3- Ji) 

bi = ^S^J2a^ 

b2 = ^S^16Js + 2J2-Ji)a^ 



(50) 



It is clear that the energy (49) is the most general expan- 
sion up to the fourth order derivatives. This is true for 
the minimal Heisenberg model (45), and this is also true 
for any extension of the Heisenberg model that does not 
include mobile fermions. Depending on relation between 
the parameters 61 and &2 the spiral is directed along a 
diagonal or along an axis of the square lattice. One can 
also tune up bi and &2 in such a way that there is a di- 
rectional degeneracy. The ground state energy gain due 
to the spin spiral is proportional to Q^. All in all this 
is very much different from the chiral scenario that we 
have considered for the "antiferromagnetic spiral". The 
effective action in the present case is 



(51) 



This is still a gradient expansion, but the theory is quali- 
tatively different from (7). The AF chiral theory (7) orig- 
inates from the t-J model and therefore the wave vector 
of the spin spiral Q is naturally small, it is proportional 
to the small doping. In the "ferromagnetic spiral" de- 
scribed by the action (51) the wave vector Q is small 
only due to accidental cancellation of different contribu- 
tions to the spin stiffness. The cancellation results in a 
very small negative stiffness. 

To be specific we consider the case 61 < 62- In this case 
Q is directed along one of the diagonals, for example 



Minimization of the energy (49) gives 



(52) 



Ps 



bi 



(53) 



B. Spin wave excitation, linear approximation 

Equation of motion generated by (51) reads 

(54) 



„. S . SE 5n 
2i—ri 

1^ 



' Sn (577^ 

This results in the following equation for n 



Dn 



(55) 



The indexes a, /3, 7 denote components of vector in the 
spin space. The differential operator in the square brack- 
ets acting on a plane wave gives a polynomial. 



D = 



a 

'J 



- -j[~Psq' + 2b,qyy + b2{qt + qt)] ■ (56) 



Looking for solution in the form 

il w (cos(Q • r + (pg),sm{Q ■ r + tpg), hq) , 



(57) 



and substituting this in Eq.(55) we obtain the following 
equations 



</5q = - [Dq - Dq] hq 
h„ = Dq - -Dq+Q - -^^q-Q 



(58) 



The fixed- frequency solution of these equations, ipq,hq cx 
exp(—iujqt), results in the following frequency 



-Dq+Q + l^Dq_Q - Dq 



[Dq - Dq] . (59) 



It is convenient to rewrite the frequency in the following 
notations 



^ofp 



(60) 



fp — \J XpYp , 
where p is the dimensionless momentum defined in (16). 



no = 


{bi^ 




Xp = 


2(P? 


+ apl) + + 2aplpl 


Yp = 




- pY + 2apjpl , 


a = 


62- 


bi 


62 + 


bi ■ 



(61) 



We remind that pi is directed along Q and p2 is perpen- 
dicular to Q. The dispersion (60) is plotted at a = 0.5 for 
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FIG. 6: (Color online) A and B: the magnon dispersion in the 
direction parallel to the spin spiral (A), p = (pi,0), and in 
the direction perpendicular to the spin spiral (B), p = (0,p2)- 
The red dashed line in Fig. B shows schematically opening of 
the gap driven by quantum fluctuations as it is discussed in 
subsection D of the present section. The "unshifted" intensity 
factors = are plotted in Figs. (C) and (D) for the 

same directions. The intensity factors at values of k where 
cjfc > O.SrJo are shown by blue dotted lines. All the plots are 
presented for the anisotropy parameter a = 0.5. 

directions parallel to the spin spiral, Fig. 6A, and perpen- 
dicular to the spin spiral, Fig. 6B. There are four points 
to note, (i) The dispersion scales as flo oc that is a 
consequence of the fourth derivatives in the action, (ii) 
The energy is zero at p = and at p = (±1,0). These 
zeroes are dictated by the Goldstone theorem, (in) The 
energy is also zero ai p ~ (0, ±1). These zeroes are not 
dictated by any exact theorem and we will see that quan- 
tum fluctuations open a gap at these points, (iv) The dis- 
persion atp ^ 1 is anisotropic, Wq « VLq [p'*/2 + 2ap\p^ . 
This is because the theory is stabilized by forth spaeial 
derivatives that "know" about the lattice structure. 

Representing Lp and h in terms of the spin wave anni- 
hilation and creation operators (cp, cp, 

h = Y^-iAqCqe-''^''*+"i'' + h.c. 
q. 

= ;^BqC,e-''^''*+"' '- + /i.c. , (62) 
q. 



with the energy transfer a; and the momentum transfer 
fc is proportional to 

l(,i^,k)=n^A\5{uj-ujk) (64) 
+ ^"^11 {^I^qK^ - ^k+q) + BI_qS{u - ujk-q)} , 

where TZ± and 7?.|| are kinematic factors that depend 
on the relative orientation of the neutron polarization 
and the plane of the spin spiral. These factors depend 
on specific experimental details and therefore we do not 
present these factors. Thus, there is an "unshifted" con- 
tribution in the neutron scattering probability with in- 
tensity proportional to Af, and two "shifted" contribu- 
tions with intensities proportional to \B^_i^q. As one ex- 
pects, the intensity is diverging at points where ujk — >■ 0. 
To characterise the intensity we plot the intensity factor 
Z^^^ ^ Yigs. 6C and Figs. 6D for directions of k 

parallel and perpendicular to the spiral. The "positively 
shifted" spectrum cok—Q is plotted in Fig. 7 for the di- 
rection of k parallel to the spiral. In the same Fig. 7 



0.5- - 2.5- 



*^ -2 ' -i ' 1 2 pj " -2 -10 1 2 pj 

FIG. 7: (Color online) A; the positively shifted magnon dis- 
persion LOh—Q for the direction of k parallel to the spiral. 
B: the intensity factor Z'j^'^ — "20 ^I-q the posi- 
tively shifted line. The intensity factor at values of fe where 
ujh-Q > O.5SI0 is shown by blue dotted line. The plots are 
presented for the anisotropy parameter a — 0.5. 

we also plot the intensity factor = "20 ^k—Q 
the shifted line. Comparing Figs. 6 C,D and Fig. 7B 
we conclude that the unshifted "line" has a significant 
amplification compared to the shifted one. 



C. Quantum and thermal fluctuations. Dependence 
of the spiral wave vector on temperature 



and applying the standard quantization procedure to the 
action (51) we find the coefficients Ap and Bp, 



'Xr. 



Bp — 



(63) 



Neutron scattering amplitude is given by the Fourier 
transform of (57) and hence the scattering probability 



Since we are interested in nonlinear corrections, Eq. 
(57) has to be written more accurately. 



n = (C^, 5^,0) COS7 -f (0, 0, 1) sin7 

= cos(<5 ■ r + Lp) 

= sin((5 - r + ip) 

h = sin 7 . 



(65) 



Equations (52) and (53) have been derived by minimiza- 
tion of the semiclassical energy. Alternatively they can 
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be obtained from the condition that the hnear in opera- 
tors term in the action (energy), 

El = -psiQxdx'P + Qydyip) 

+ 2blQxQy [QydxSP + Qxdyif) 

+ 2b2 {Qld^ip + Qldy^) , (66) 

is zero, Ei = 0. Using (65) and expanding Eq. (49) 
up to cubic terms in perturbations h and tp we find the 
foUowing cubic terms 

(Vn)2 ^ -2h'^{Q-V)ip (67) 

{d^^dyfif 2{Q^,d^ip + Q^d^,'f)[{dfj,if){d„(p) 
- Q^,Q,h^ + 2d^ihd,h)] . 

There is no summation over fi and v in (67), so one can 
set fi = X, i> = y or = X, V = X, etc. To calculate the 
tadpole, Fig. 4, with interaction (49), (67) we note that 
the following expectation values are nonzero 

M = {{d^MidxV)) = {{dy^){dyV)) 



2AI (n. ■ ' 



Here 



(68) 



(69) 



is the usual Bose distribution. We remind also that the 
dimensionless p is related to q by Eq.(16), and note that 
{d^{n;,d„n,_)) = 0. 

When calculating the tadpole the Hartree-Fock decou- 
pling in Eqs.(49) and (67) gives two kind of terms. Terms 
of the first kind proportional to (n^) lead to the following 
renormalization 



Ps P3(l-(n2)) 
fei ^ h,{l- {nl)) 
62 62(1 -(n^)) 



(70) 



Obviously, this renormalization does not change Q given 
by Eq.(53). Terms of the second kind are proportional to 
M and N defined in (68). They generate the following 
fluctuation correction (tadpole) to the effective energy 

= 2M{hi + Zh2){Qxdx'4^ + Qydy^) 

+ 4J\fbi{Q^dy(p + Qyda^ifi) . (71) 

The renormalized spiral wave vector Q has to be found 
from the condition E1+E3 = 0, see Eqs. (66), (71). This 



gives the following correction to Q 
f = -^E^^ [3,? + (l + 2a)9|] 




2a)p\ 



(72) 



We remind again that the momentum p in this equation 
is dimensionless, p = q/Q. The quantum fluctuation part 
of SQ, that comes from 1/2 in (rip -|~ 1/2), is ultraviolet 
divergent as (in 3D it would he p^). Therefore it must 
be renormalized to zero. The temperature dependence of 
Q is given by the Tip-term that is left after the renormal- 
ization. 

To understand better the quantum renormalization 
procedure we rewrite (66) and (71) having in mind that 
Qx = Qy = Ql V2, 

El = [-psQ + {ai+a2)Q^]diip 

E3 = [2Miai+3a2)Q + 4AfaiQ]diip . (73) 

Using (68) one can easily calculate A^, A/" at T=0, 



M = 



N 



a 

2S 
2S 



167r 



A 



A 



(74) 



where A ^ 1/a is the ultraviolet cutoff and 
/2, /4, /ei 52, 547 56 are some functions of the anisotropy 
a that can be easily calculated. Comparing (74) with 
(73) we see that the A^ term renormalizes ps and the 

A^Q^ terms renormalize (ai +02). The Q^'^ln^^^^ and 

Q'^ terms give contributions to the energy density of the 

order of ~ ^{Qa^Q'^. So, the In (^^^ and Q** terms 

in (74) have to be neglected since in (49) we keep only the 
terms ^ JQ^ and ~ J{Qa)'^Q^ neglecting higher deriva- 
tives. 




FIG. 8: Temperature dependence of the integral Ft in Eq.(75) 
for Q{T). The value of a is 0.5 

Temperature dependence of the spiral wave vector is 
determined by the "rip-part" of Eq. (72), 



SQ 



Ft 




(75) 



[3pl + (1 + 2a)pl] Up 
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The dimensionless integral Ft is plotted in Fig. 8 versus 
T/fig for oi = 0.5. Temperature reduces the value of Q 
and this is a sizable effect. 



D. Quantum correction to the spin-wave 
dispersion. Opening of the gap at q = (0, ±Q) 

There is the tree leg vertex in the theory shown in Fig. 9 
The vertex is determined by Eqs.(49) and (67). There 




FIG. 9: The tree leg vertex 

are three contributions to the vertex. The first contribu- 
tion comes from the 7 ^ip terms shown in Fig. 10 left. 
The second contribution comes from the ip — >■ 77 terms, 






FIG. 10: Three different contributions to the 3-leg vertex 

Fig. 10 middle, and the third contribution comes from the 
ip — > ipip terms, Fig. 10 right. A straightforward calcula- 
tion gives the following answers for these contributions 



F 

Fi 

F2 
F3 



= F1+F2+F3 

= 2iQ{ai + 02) X {pi [pI+pI{1 

+ h [ll + ll{l + 2a)] AgBiAp} 

= 2iQ{ai + 02) X qi [ql + q^^l 4 

= -2iQ(ai 4- 02) X {Sqipih 

+ (1 + 2a) [qiP2h + Pi92^2 + ^1^292]} BqBpBi 



(76) 

+ 2a)] AqBpAi 
2a)] BqApAi 



Importantly, for an on s/ie// scattering process, q ~ p + l, 
Ljq ~ Ljp + oji, the vertex satisfies the Adler's relation, it 
vanishes if one of the momenta equals to or Q. 

To understand the structure of the quantum correc- 



tion to the spin-wave dispersion we calculate the 7 
polarization operator shown in Fig. 11 
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FIG. 11: The 7 — 7 polarization operator 



2TlTl{LL>q/2-k 



'•^q/2+kY 



(77) 



The polarization operator can be represented in the fol- 
lowing form 



P2-'{uj) = -A 



Q^(ai 



0,2, 



S 



c? X integral 



(78) 



The integral is quadratically ultraviolet diverging. Simi- 
larly to (74) only ultraviolet diverging contributions are 
important with the accepted accuracy. Calculation of di- 
verging terms in the integral is very simple. For example 
the result at a = 0.5 is 



integral 
x(8.7Q2g2 



(5.2g2 + i.5g|)^ 
4.7qt - 



5.0Q2q2 



ln(A/Q) 

27r 
1.492' 



(79) 



lOqfqi) 



The polarization operator determines the 7j7q term in 
the effective action. The coefficient A"^ in (78) gives nor- 
malization of the 7-ficld and powers of momentum in (79) 
should be replaced by corresponding derivatives of this 
field. All in ah Eq. (78) together with Eq. (79) resuh in 
the following effective energy 



0.2} 



5.2 



(Aa)2 



47r 



.7(Qa) 



2 ln(A/Q) 
27r 



(80) 
(917)' 



1.5 



[Kaf 



4tt 

ln(A/Q) 



5.0(Qa) 



2ln(A/Q) 



2tt 



{d2lf 



. 2^ a^[4.7{dhf + lA{dh? 

+ 10(^1 927)']} 

All momenta in the physically interesting region arc of 
the order of Q. Hence, the typical value of the bare en- 
ergy density (49) is E'bare ^ The logarithmic terms 
in (80) are of the order of ~ ^(Qa)^, so they must be 
neglected within the accepted long-wave-length approxi- 
mation, Qa <^ 1 Only the power divergent terms in (80) 
are significant within the long-wave-length approxima- 
tion, hence 



Q2(ai+a2) (Aa)2 



The structure of SE, 



is the same 



5.2(9i7)2 + 1.5(927)'] (81) 



The most impor- 
tant issue is that the polarization operator correction (81) 
contains different coefficients in front of the derivatives 
in different directions. To satisfy the Goldstone theorem 
we must introduce a contra-term —^6ps{'Vn)'^ to renor- 
malize out the longitudinal derivative, (^17)^. Because of 
the mentioned different coefficients the term (927)' is not 
zero after the renormalization. Therefore, we necessarily 
open a gap in the dispersion at the points q = (0,±Q). 
The gap comes from fluctuations at the ultraviolet cut- 
off, therefore the effective field theory docs not allow to 
calculate the gap. An estimate for the gap follows from 
(81), 



A- 



ttS 



(82) 
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Here we take into account that A ^ 1/a. The gap open- 
ing is shown schematically in Fig. 6B. Due to quantum 
fluctuations the gapless red dashed part of the disper- 
sion disappears and the dispersion becomes gaped. It is 
worth noting that the mechanism of the gap opening due 
to quantum fluctuations is similar to that in the spin- 
stripe phase of the Ji — J2 model.^*'^^ 



E. Discussion and conclusion of the "ferromagnetic 
spiral" section 

We have formulated the effective field theory describ- 
ing a spin spiral close to a coUincar ferromagnetic state, 
i.e. the wave vector of the spiral is small compared to in- 
verse lattice spacing, Q <^ n/a. The problem is inspired 
by recent data^ on FeSrOa and FeCaOs. We assume 
that there are only localized spins, no mobile spin carri- 
ers. This assumption results in the effective action (51) 
that contains forth spacial derivatives. Using the effec- 
tive action we calculate the spin- wave excitation spectra, 
see Figs. 6 A,B. The spectra have three zero modes at 
q = 0, ztQ dictated by the Goldstone theorem. In the 
semiclassical approximation there are also two additional 
zero modes that are not dictated by any exact theorem. 
We demonstrate that quantum fiuctuations open gaps at 
these points. 

We have also calculated the temperature dependence 
of the spiral wave vector, Q = Q{T), see Eq.(75) and 
Fig. 8. 

The typical energy scale for the spin-wave excitation 
with momentum q ^ Q is Hq (x Q"^, see Eq.(61). The 
dependence oc is a consequence of the forth order 
derivatives in (51). Based on data^ we estimate that 
flo ~ lOmeV in FeSrOa and FeCaOs. Therefore an ex- 
ternal magnetic field of the order of several Tesla can 
substantially modify the spin spiral making it conical. 

Last but not least. The intensity of inelastic neutron 
scattering grows monotonously at large q, see Figs. 6 B,C. 
Within the accepted model we do not see a way to sup- 
press intensity outside of the region q ^ Q- We did cal- 
culate broadening of magnons due to the decay 1—^2, 
see Fig. 9, and found that the broadening is very small. 
Introduction of mobile spins in addition to the localized 
ones can change the situation. In this case a redistribu- 
tion of the quasiparticle weight between two brunches of 
the dispersion similar to that for the AF spin spiral is 
possible. However, discussion of this effect is outside of 
the scope of the present paper. 



V. CONCLUSION 

In the present paper we derive and analyze effective 
quantum field theories for helical spin magnets close to 
antiferromagnets and to ferromagnets. We assume that 
the wave vector of the spin spiral is small compared to the 
wave vector of the lattice, Q ^ 27r/a. The noncoUinear 



structure is caused by spin frustration. The field theory 
describing the "antiferromagnetic spin spiral" is inspired 
by physics of cuprate superconductors. In this case the 
wave vector of the spin spiral is naturally small because 
of the small doping x. The field theory describing the 
"ferromagnetic spin spiral" is inspired by the recently 
studied FeSrOa and FcCaOs compounds. In this case the 
wave vector of the spin spiral is small due to accidental 
fine tuning of parameters. 

The effective action describing the AF spin spiral (7) 
contains two bosonic fields, n corresponding to the stag- 
gered spin, and ^ corresponding to pseudospin. The ac- 
tion is quadratic in spacial derivatives, contains the sec- 
ond order time derivative of n and the first order time 
derivative of the spinor field describing ^ in the CP"'^ rep- 
resentation. The action is space and time covariant, the 
covariance is a consequence of the gauge invariance of 
the t-J model which is the origin of the effective action. 
The "hourglass" dispersion of magnons, see Figs. 2, 3, 
is a direct consequence of the action. The two branches 
of the "hourglass" dispersion originate from two differ- 
ent degrees of freedom, spin and pseudospin. The lower 
branch of the "hourglass" has strongly anisotropic mo- 
mentum dependence with two zero modes at the Gold- 
stone points q — ±Q, while the upper branch is prac- 
tically isotropic. This behavior agrees with the typical 
dispersion observed in cuprates.^^^^^ Dependent on pa- 
rameters of the action the AF spin spiral can be either 
static or dynamic. There is a Quantum Critical Point 
that separates the static and the dynamic regimes. The 
spin spiral itself is characterised by the pseudoscalar chi- 
ral order parameter independently whether this is static 
or dynamic case. The chiral order parameter changes 
sign under time inversion. 

The effective action of the ferromagnetic spin spiral, 
Eqs. (51) and (49), contains only one bosonic field n de- 
scribing spins. The action contains the second and the 
fourth order spacial derivatives as well as the first order 
time derivative of the spinor field describing n in the CP^ 
representation. Dependent on parameters of the action 
the wave vector of the spin spiral Q can be directed ei- 
ther along a diagonal or along an axis of a square (cubic) 
lattice. The effective action results in the spin-wave ex- 
citation spectra shown in Figs. 6 A,B. The spectra have 
three zero modes at q = 0, ±Q dictated by the Gold- 
stone theorem. In the semiclassical approximation there 
are also two additional zero modes. Quantum fluctu- 
ations open gaps at these points. The wave vector of 
the spiral has a sizable temperature dependence that has 
been calculated. 
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